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The strong coupling limit (j3 gaug e = 0) of QCD offers a number of remarkable research possi- 
bilities, of course at the price of large lattice artifacts. Here, we determine the complete phase 
diagram as a function of temperature T and baryon chemical potential jig, for one flavor of stag- 
gered fermions in the chiral limit, with emphasis on the determination of a tricritical point and 
on the T m transition to nuclear matter. The latter is known to happen for jUg substantially 
below the baryon mass, indicating strong nuclear interactions in QCD at infinite gauge coupling. 
This leads us to studying the properties of nuclear matter from first principles. We determine 
the nucleon-nucleon potential in the strong coupling limit, as well as masses toa of nuclei as a 
function of their atomic number A. Finally, we clarify the origin of nuclear interactions at strong 
coupling, which turns out to be a steric effect. 
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1. Model and motivation 

We study lattice QCD with one species of staggered fermions at infinite gauge coupling [j|]. 
The partition function 

Z(m q ,ii) = J @U@X@Z e 5p , (1.1) 

depends on the staggered quark mass m q and the baryon chemical potential jXg- Since j3g aU ge = 0, 
we only have the fermionic action Sp with 

Sf= £ VxfXx U x jXx+v-U^vX x -v +2m q ^XxXx (1.2) 

x,v=l,4 x 

where r\ x % = (— 1)Ep<v*p 5 r\ x j = 1. The quark chemical potential ji. and an anisotropy 7 are intro- 
duced by multiplying the time-like gauge links U ±% by 7exp(±a/x/7) in the forward and back- 
ward directions, respectively. The anisotropy allows to vary the temperature continuously, via 
T=aT l 'f/N % at infinite coupling 

When m q = 0, our model is invariant under global U(1)a x U (l)g transformations 

X(x) -> / £ W^ +es) z(x), £(x) -> ^(i)/ e W e "- e »)Vi, (1.3) 

where e(jc) = (— l)^* Xv . While £/(l)g is responsible for baryon number conservation and remains 
unbroken, U(1)a represents the chiral symmetry of the model and can break spontaneously. 



Owing to the absence of the gauge action in Eq. (1.2), the link integration factorizes in Eq. (1.1) 



and can be done analytically (4)]. The degrees of freedom are now mesons and baryons. Carrying 



out the Grassmann integration, the partition function Eq. fll.ip becomes for m q = that of a dimer- 
loop model |jl|] : 

{„ ,.<■} v.v n *t- c 

with the constraint that mesonic links with occupation number n x $ = 0, ..,3 attached to a site x 
satisfy Y,±u. n x,±fi = 3- Alternatively, a site can be traversed by a self-avoiding baryon loop C. 
The weight of such a loop C is given by w(C) = p(C)^ Ni ^ exp (ik\iaN x /Y). Here, N%(C) is the 
number of links on the loop in the time direction, k is its winding number in this direction and 
p (C) = ± 1 is a geometry-dependent sign, which can be negative even when pL = 0. Karsch and 
Mutter [p removed the sign problem present already at ji. = by analytically resumming pairs of 
configurations. When /I > the remaining sign problem is mild (see next Section), allowing us to 
simulate large enough lattices that the phase diagram of the model can be determined reliably. 

There exists a plethora of mean-field results dating back to the early days of lattice QCD []|] 
and continuing up to now [^, ^J, with the inclusion of NLO and NNLO corrections to the infinite 
coupling limit [Q], and statements about a quarkyonic phase [^]. These approximate predictions 
should be checked against numerical simulations using an exact algorithm. 

Moreover, lattice QCD should provide the means to study nuclear physics from first principles. 
In lattice QCD at weak coupling, properties of single hadrons [fill], nuclear scattering lengths and 



potentials [11, 12, 13] or two-and three baryon systems [14] are being studied. Going to higher 



nuclear matter density is mostly hindered by a severe sign problem. Contrary to that, the strong 
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Figure 1: Finite-size scaling of the chiral susceptibility near the phase transition, with (a) d — 3, 0(2) 
exponents at a\i = and afi = 0.3, and (£>,top) tricritical (mean-field) exponents at ajj. = 0.33. Lines are 
drawn to guide the eye. For a/x = 0.36 (b, bottom), we instead rescale the y-axis with L-7l v = L- 3 and 
observe that Xcr/L 3 goes to a constant for T < T c (fi). The solid lines are fits derived from the first-order 
ansatz Eq. (2.4). All results have been obtained for N T = 4. Errors are of the size of the symbols. 



coupling limit offers not only the possibility to study nuclear matter at higher density thanks to the 
mild sign problem, but also provides an intuitive understanding, owing to the point-like nature of 
baryons and pions in this limit. The absence of earlier numerical studies of the strong coupling 
version of nuclear physics is mostly due to algorithmic issues present in the early studies []IJ [l5|] . 
Decisive progress occurred with the introduction of the worm algorithm JTq], which has been 
adapted for strong coupling SU (2) and U (3) lattice theories in [17], enabling efficient Monte Carlo 
sampling even in the chiral limit m q ^0. We extend this approach to SU (3). 



2. Phase diagram 



Since we consider here the case of a massless quark, the chiral symmetry U (\)a Eq. (1.3) is 
exact but spontaneously broken at small (T,fi), with order parameter When pL = 0, a mean- 

field analysis predicts symmetry restoration at aT c = 5/3, whereas the Monte Carlo study of [Q] on 
N t =4 lattices, extrapolated to m q = 0, finds aT c = 1.41(3). In order to determine T c we performed 
numerical simulations at a\i = 0,am q = 0. Our main observable is the chiral susceptibility 



1 92 i v 

V ami 



(£wWw(o)), (2-1) 



where we used the fact that in a finite volume the order parameter (iffy) vanishes identically. 
Correspondingly, % a will not show a peak at the transition but finite size scaling (FSS) still applies, 
so that for a system of size L 3 x N x at a "reduced temperature" t = I — T /T c {or \ — \ij \L C ): 

X a ~iyl v x{tl} ,V )=L y ' V (ao + 0{tL l ' v )) . (2.2) 
In principle, this allows us to determine the exponents y and v, which should be those of the 



d = 3, 0(2) universality class [2C]. In practice however, the numerical values are close to those 
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(a) (b) 

Figure 2: (a) (n,T) phase diagram of 1 -flavor strong coupling QCD with massless staggered fermions 
(N T = 4). (b) a 4 Af = — y 2 log ((sign) )/(L?N T ) versus a\i at T — Tjcp- Lines are drawn to guide the eye. 



of the Z(2) universality class and we point to [^I|] for a more advanced treatment in the case of 
U(3). Here, we simply assume d = 3, 0(2) exponents and show that Xa collapses on a universal 



curve when rescaled according to Eq. ( ]2.2| ). For example, Fig. |l(a)| (top) shows the collapse for 
/i = 0,N r = 4 and several values of L, with aT c fixed to 1.402 which is the crossing point of % a jTjl v 
when plotted as a function of T. Using this strategy, we find aT c = 1.319(2), 1.402(3), 1.417(3), 
respectively, for N t = 2,4,6, indicating an Af T ->oo limit about 15% smaller than the mean-field 
prediction. 

For jU > the model has a sign problem. We have measured the "average sign" 

(sign) = |- = exp (-AfV/T) = exp(-a 4 A/L 3 Ar T / y 2 ) , (2.3) 
Z ll 

where Zi corresponds to the partition function Z of Eq. (|L4]), but taking the absolute value of the 
weights. The free energy density difference (a 4 Af) is a measure of the severity of the sign problem. 
Fig. [2(b)| shows (a 4 Af) as a function of a\i, for several L 3 x 4 lattices at T = 0.931 a^ 1 7tcp (see 
below), using the analytic resummation prescription of Karsch and Mutter [|l|] which removes the 
sign problem at jU = 0. Af is nicely volume-independent, vanishes at /I = 0, starts <=c p} (see inset), 
and peaks slightly past the phase transition. Note the very small magnitude ^(10 -4 ) - compared 
to (a 4 Af) ~ &(\) expected when using the standard approach of integrating over the fermions first 
- which allows us to simulate 16 3 x 4 lattices with (sign) > 0.1. 

For the available volumes, we may then follow the critical line as a\i increases, monitoring 
the collapse of % a using the appropriate critical exponents. Expectations are that the second order 



0(2) transition will turn first order at a tricritical point (TCP) for some nonzero /x. From Fig. |l(a) 



(bottom), we see that for ajl = 0.3, % still obeys 0(2) scaling behavior. With a slight increase 



to a\i = 0.33 Fig. |l(b)| (top), a satisfactory collapse requires tricritial (mean-field) exponents y t = 
1, v t = 1/2. Under a further small increase to a\i = 0.36 no such finite-size scaling collapse can be 
achieved. Instead, the transition is first-order: (i) the distribution of the baryon density shows two 
peaks, whose areas become equal at aT c = 0.844(1); (ii) for T < T c , Xa/L 3 becomes L-independent 



(see Fig. 1(b) (bottom)); (Hi) on either side of T c , % c is well described by differentiating the two- 



4 



Nuclear physics at strong coupling 



phase Borgs-Kotecky ansatz [p3[] 



z(T) = eX p(-^/ 1 (r))+ ceX p(-^/ 2 (r)) (2.4) 

with f\z{T) = fo ± Cf (r — r c ), from which the solid curves can be obtained. Thus, from the avail- 
able data we conservatively conclude that (a^TCP) a ^TCp) = (0.33(3), 0.94(7)). This should be 
compared with the analytic prediction (0.577,0.866) of [0|. The rather large difference in un- 
derlines the ff (l/d) accuracy of a mean-field treatment, and justifies a posteriori our Monte Carlo 
study. In Fig. [2(a)| we summarize our findings for the phase diagram. In spite of its resemblance 
to the expected deconfinement transition in massless Nf = 2 QCD, here the two phases are both 
confining, with point-like mesons and baryons, and so the phase transition is to dense, chirally 
symmetric nuclear matter. At T = the baryon density j umps from to 1 , a saturation value caused 
by the self-avoiding nature of the baryon loops, which itself originates from their fermion content. 
Using the baryon mass to fix the lattice spacing, this represents about 4 nucleons per fm 3 , around 
25 times the real- world value. An intriguing feature of this T = transition - and an important mo- 
tivation for this study - is the value of /x| ntlcal , which both mean-field [j|] and an early Monte Carlo 
study [jl]] find much smaller than the naive threshold value rag. However, the ergodicity of the sim- 
ulations of [p]] was questioned in Ji3|], which was found to be justified in [24]. This motivated us to 



redetermine /ij nt (r = 0) using an improved method inspired by the "snake" algorithm [25]: When 
two phases coexist, the free energy AF /T necessary to increase by a "slice" Lx Lx a the volume 
occupied by dense nuclear matter can be decomposed into L? elementary contributions, looking 



generically like Fig. 3(a), where one additional static baryon is attached to 3 neighbors. We mea- 
sured the free energy AF /T of this elementary increment on a large 8 3 x 16 lattice, and obtained 
aAF = a/J.g nt = 1.78(1), rather close to both mean-field predictions [||] and Monte Carlo extrapo- 
lations This we compare to artis which corresponds to the difference in free energy measured 
at jJ. = 0,T ps by extending a static baryon world line. We find am% = 2.88(1), consistent with 
HMC simulations [18], and again in agreement with mean-field [|5|] and large-A^ [jl^] predictions 



but indeed much larger than a\i\ ,. As already recognized in [|26|], the reason that /i™ < ra# must 



then be the presence of a strong nuclear attraction, to which we now turn. 

3. Nuclear Physics at infinite gauge coupling 

Since our baryons are point-like, there is no conceptual difficulty in defining the nuclear po- 



tential Vnn(^), unlike in the real world We measured Vnn using again the "snake" algorithm, 
this time extending little by little in Euclidean time the worldline of a second baryon at distance 
R from the first. The result is shown in Fig. 3(c)| . Aside from the hard-core repulsion, there is 



indeed a strong nearest-neighbor attraction, a slight repulsion at distance ay/l, and almost no in- 
teraction beyond that distance. Vnn is similar qualitatively to what is expected in the real world, 
with competition between attractive a exchange and repulsive ft) exchange. The depth of the min- 
imum ~ 120 MeV and the corresponding distance ~0.6 fm are quantitatively plausible fl28|]. This 
nearest-neighbor attraction also explains a posteriori the value of /xj nt : each baryon added to the 
dense phase binds with 3 nearest neighbors, which reduces the increase in free energy from ara# to 
only a(ra# + 3Vnn(o)) ~ 1-7, consistent with a/i| nt - 
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Figure 3: (a) Adding a baryon to grow an additional layer of bulk nuclear matter: each new baryon binds 
to 3 nearest-neighbors, (b) Building a first layer of nuclear matter inside the hadron gas, thus creating two 
interfaces: each new baryon binds to 2 nearest-neighbors, (c) Energy of a second static baryon at distance 
R from the first, i.e., (mb + Vnn(2?)), where Vnn(^) is the nuclear interaction potential. The horizontal band 
indicates the mass of an isolated baryon and corresponds to Vnn =0. At R = the potential is infinitely 
repulsive. 

Similarly, we can predict the T = surface tension of nuclear matter: in a periodic cubic 
box, when building a first "slice" of nuclear matter with two interfaces in the dilute phase, each 
new baryon binds with only 2 nearest-neighbors (Fig. 3(b) ) instead of 3 in the bulk (Fig. 3(a)| ), 
thus increasing its free energy by |Vnn(a , )| for an increase of 2a 2 in the interface area, yielding 
(7~^|VnnO)|- 

This large interface tension has an impact on the stability of nuclei of various sizes and shapes: 
for a given atomic number A, those with a shape close to a sphere (or a cube) will have a smaller 
mass. Using the same variant of the "snake" algorithm, we have added baryons, one by one, to form 
such nuclei while measuring the successive increments in free energy. For A = 2 our "deuteron" 
binding energy is about 120 MeV: the real-world binding energy of ~2 MeV results from delicate 
cancellations which do not occur in our 1 -flavor model, and the binding energy remains of the same 
magnitude as the depth of Vnn- For larger A, the resulting Fig. |4(a) does indeed show increased 
stability for nuclei having square (A = 4), cubic (A = 8) or parallelepipedic (A = 12) shapes. Other 



"isomers" with different shapes, studied exhaustively for A = 4 and sketched Fig. |4(b)| , have clearly 
larger masses. Moreover, the average mass per nucleon is well described by the first two (bulk and 
surface tension) terms of the Weizsacker phenomenological formula: 



m(A)/A = + (36n) i/3 a 2 aA 



-1/3 



(3.1) 



where a is set equal to Vnn (a) I in the Figure. The next higher-order terms in this formula come 
from isospin and Coulomb forces, which are both absent in our model. 



4. Summary and Remarks 

An interesting aspect of our study is the origin of the nuclear interaction. The nucleons are 
point-like and self-avoiding, so that only the hard-core repulsion is explicit. There is no pion 
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Figure 4: (a) Mass per nucleon of A — 1, .., 12 nuclei. For A — 3,4 all possible geometric isomers are 
included. The solid line shows the parameter-free Bethe-Weizsacker Eq. (3.1), with the surface tension a 

_2 

set to ^|VNN( a )|- (b) Corresponding nuclear geometries in order of increasing mass, (c) Energy density of 
the pion cloud as a function of the Euclidean distance to a static baryon. 



exchange. In a way reminiscent of the Casimir effect between two neutral plates, the interaction 
proceeds by the rearrangement of the pion bath caused by the excluded volume of the nucleon. 



This rearrangement is visible Fig. |4(c)| for one nucleon: at a neighboring site, the three pion lines 
attached to each site have fewer options and orient more often along the Euclidean time, which 
increases the pion energy. In fact, the nucleon mass awi£«2.88 can be decomposed into a bare 
mass 3 — 3/4 = 2.25, which is the energy increase "inside" the nucleon and can be assigned to 
the three valence quarks, and an energy increase « 0.63 in the surrounding pion "cloud". When 
two nucleons are next to each other, the latter increase is limited to 10 nearest-neighbors instead 
of 2 x 6, which explains the attraction between them (in sign and roughly in magnitude). This 
excluded volume or "steric" effect is thus the origin of the nuclear potential, and ultimately of 



nuclear stability, in our model. In real QCD, the pion density is not constrained as in Eq. (|1.4|). 
Nevertheless, it is going to be high at temperatures T ^m n [ |29| ] and one should expect the same 
steric effect to enhance nuclear attraction at such temperatures. 

To summarize, in a crude model of QCD, 1 -flavor massless lattice staggered fermions at strong 
coupling j3 gaug e = 0, we have been able to obtain the complete phase diagram and derive the strong 
coupling version of nuclear interactions and nuclear masses from first principles, uncovering a 
simple, but universal, steric origin of the nuclear interaction. This model can be improved in many 
ways. One simple modification consists of giving a non-zero mass to the quarks: the nuclear inter- 
action will weaken as the pion mass is increased, in a way which can be compared with effective 
field theories. Less simple but feasible improvements include introducing isospin with a second 
quark flavor, and measuring the ^(j6) correction as done analytically in [|(], 31, §§]. These will 
bring our model much closer to real QCD. 
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